Trap and population imbalanced two-component Fermi gas in the BEC limit 
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We study equal mass population imbalanced two-component atomic Fermi gas with unequal trap 
frequencies [ijJ^^ ^ uji) aX zero temperature using the local density approximation (LDA). We consider 
the strongly attracting Bose-Einstein condensation (BEC) limit where polarized (gapless) superfluid 
is stable. The system exhibits shell structure: unpolarized SF— >polarized SF^normal N. Compared 
to trap symmetric case, when the majority component is tightly confined the gapless superfiuid 
shell grows in size leading to reduced threshold polarization to form polarized (gapless) superfluid 
core. In contrast, when the minority component is tightly conflned, we flnd that the superfluid 
phase is dominated by unpolarized superfluid phase with gapless phase forming a narrow shell. The 
shell radii for various phases as a function of polarization at different values of trap asymmetry are 
presented and the features are explained using the phase diagram. 

PACS numbers: 03.75.Ss, 03.75.Hh, 67.85. Lm 
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INTRODUCTION 

Ultracold atomic Fermi gases present a unique oppor- 
tunity to study the exotic pairing phases where the effec- 
tive interaction is tunable via Feshbach resonance and 
population of each spin component can be controlled. 
These studies with two-component Fermi gas include the 
Bardeen-Cooper-SchriefFer (BCS) to Bose-Einstein con- 
densation (BEC) crossover P, Q with equal population 
two-component Fermi gas and the effects of population 
imbalance on the superfluid state [2, iS, H, S H B Q • 
For Fermi gas with population imbalance, various pairing 
scenarios arc proposed; Fulde-Ferrel-Larkin-Ovchinnikov 
phase FFLO 10, 11], breached pairing 12], phase sepa- 
ration [l^ and pairing with deformed fermi surface [3]. 

The present studies of equal mass superfluid Fermi sys- 
tems involve same trapping potentials for the two com- 
ponent. The zero jl^ and finite temperature [IB] phase 
diagrams of the equal mass population imbalanced sys- 
tem taking into account various pairing phases with im- 
plications to experiments are well understood. The trap 
imbalanced is naturally realized with Fermi mixture with 
unequal masses where each component experiences differ- 
ent potential due to the mass difference. The ground 
state properties for this system have been studied in 
Ref. 
Ref. 
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18, 3, 23]- However, it was recently proposed in 
that even equal mass Fermi mixture can admit 
trap imbalance and the system was studied with popula- 
tion balance [2^ and small trap imbalance [2^. 

Motivated by recent experiment performed in the BEC 
regime of interaction exploring Bose- Fermi mixture 0, 
we consider the trap and population imbalanced (tj| ^ 
ojl) Fermi mixture in this regime. The phase diagram 
in this regime becomes richer with the existence of the 
gapless (polarized) superfluid also referred to as breached 
pair phase with one fermi surface. We study the effects 
of trap asymmetry on the shell structure as function of 
polarization. The shell structure, in general, consists of 



three phases: unpolarized superfluid (BCS SF) at the 
center, the breached pair phase with one fermi surface 
(BPl) forms the intermediate shell finally surrounded 
by polarized normal (N). We do not consider the FFLO 
phase [13, [nl where cooper pail carries finite center-of- 
mass momentum. This phase is stable within very nar- 
row window of the applied chemical potential difference 
in the BCS regime. 



FORMALISM 

The system we consider is a trapped cloud of two- 
component Fermionic mixture confined by harmonic 
isotropic potential Vra {f) where r measures the distance 
from the trap center. The Fermi atoms have unequal 
population of the two pseudo-spin (hyperfine) states and 
interact via point-contact s-wave interaction. To pro- 
ceed further we start with system without the trap and 
later include it under local density approximation. The 
Hamiltonian density {h = \) for the fcrmions in this case 
is given by 



_ fi^) «',(r)+g^'|(r)^'|(r)«'i(r)*t(r), 

(1) 

where ^'^.(r) creates a pseudospin-a fermion at position 
r, £ka — fc^/2mcr; and mg. are the chemical potential 
and mass for pseudospin-cr component with cr =|, |. g 
is the bare coupling constant characterizing interparticle 
interaction. It is related to s-wave scattering length a of 
the system by the Lippmann-Schwinger equation 
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where the average kinetic energy £k ~ (sik + £2k)/2, V 
is volume and reduced mass rh = m|m|/(TOt + TO|). The 
pairing Hamiltonian can be diagonalized by Bogoliubov 
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transformations using thermofield dynamics techniques 
as described in [2l|, [l^l • This formahsm also accounts for 
the polarized superfluidity when two components have 
unequal population or masses. This leads to thermody- 
namic potential density, 
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where we have defined gap or order parameter A(r) — 
—g (\I'|(r)4'|(r)) . Introducing the chemical potential fi = 
{nl + /2 and the chemical potential difference or the 
Zeeman field h = — ii^) /2, we define quasiparticle 
energies E„ ^ E±6^ with = ek - -E^k = + A^, 
for equal mass case 5^ = —h and (3 = 1/T. 

Note that the there are two branches for the quasipar- 
ticle energies and lead to gapless modes when Ek^ — 
assuming |-fermions to be majority. At zero temperature 
the excess fermions are accommodated in negative quasi- 
particle energy states. This gapless phase is referred to 
as breached pair phase with one Fermi surface (BP I) as 
it is stable only for /it < and hence possesses one Fermi 
surface. In the deep BEC regime this phase can be un- 
derstood as mixture of composite bosons and fermion 
quasiparticles [lit . 

The gap equation is given by the condition of ex- 
tremum of thermodynamic potential density ^ = 0. 
The average number density and density difference equa- 
tion respectively are given by n = — and m ~ 

Before considering the trapped system, we construct 
the zero t emp erature phase diagram in grand canonical 
ensemble 2J| of fixed n and h. 

We start with h = and find the point where super- 
fluid state make continuous transition to vacuum state 
of molecules. This value of ^ is denoted by Hc- For small 
h < hm this behavior i.e. superfluid-to- Vacuum persists 
and leads to vertical phase boundary in the diagram. 

Next, we start with < and increasing h. Here sys- 
tem evolve from vacuum state to polarized normal state 
as /i| = II + h is now positive quantity leading to finite 
population of the "f-fermions. There cannot be superfluid 
phase here as /i < /ic- This leads to Vaccum-to-Polarized 
N phase boundary. 

Similarly we start with fi > fic with increasing h. 
Here system makes continuous transition to Breached 
pair state as superfluid starts to admit finite polariza- 
tion. The superfluid-to BPl boundray is calculated by 
numerical comparison of the thermodynamic potential 
in the respective states. As we further increase the h, 
the BPl eventually make transition to polarized normal 
state. However, depending on the value of fj, the BPl-to- 
Polarized transition can be first or second order. The tri- 
critical point where first and second order transition meet 
is indicated in the phase diagram. Also the BPl-Normal 
first order curve intersect the superfluid-BPl curve at 
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FIG. 1; (Color online) The zero temperature phase diagram 
for (/cFffl) 1 = 2.0 showing unpolarized superfluid (BCS SF), 
polarized superfluid (BPl), vaccum and polarized normal (N) 
phases. The upper (blue) dot denotes the point beyond which 
BPl state ceases to exist. The lower (red) dot represents the 
tricritical point. The dashed (red) line indicates the second- 
order transition between unpolarized Sf and BPl phase. The 
dot dashed (black) line indicates flrst-order transition between 
SF to polarized N state and BPl to normal above and below 
the upper (blue) point respectively. 



large fj,. Beyond this intersection point, BPl ceases to 
exist and there is direct first order superfluid-to-Normal 
transition. 



TRAPPED FERMI MIXTURE 

We next consider the trapped fermions confined by 
harmonic isotropic potential. This can be taken into 
account via local density approximation (LDA) where 
trapped system is treated as locally uniform with lo- 
cal chemical potential Ha-ir) — — Vto-- Mo- is actual 
Lagrange multiplier constraining number of atoms and 
Vtct — ^m^Q.^r^ with u]„ as trapping frequency for com- 
ponent a. Furthermore, the chemical potentials of each 
component at a given point in the trap can be written as 



/ii(r) = n{r) + h{r) 
H2{r) = li{r) - h{r) 



(4) 
(5) 



where ^(r) = /i — Vt(?') and h{r) — h — SVrir) where 
VT[r) = {VT] + VTi)/2 and 5VT{r) = (Ftt - Vtj,)/2. The 
quantities /i and h are determined by enforcing particle 
number constraints namely total number of atoms and 
population difference respectively. 

In terms of n{r) = n]{r) +n|(r) and m(r) = n]{r) — 
n|(r), the total number of atoms TV and population im- 
balance AA^ are given by A^ = / (Prn{r) and AA^ = 
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FIG. 2: (Color online) (a) The three radii Rfi (outer bound- 
ary of unpolarized superfluid), Rtf (outer boundary of BP-1 
phase) and Rf2 (outer boundary of N phase) plotted as a func- 
tion of polarization P at trap asymmetry parameter r] = and 
(kpa)^^ = 2.0. (b) The molecular density n™, and magnetiza- 
tion m plotted against radius r measured in units of kp and 
Rtfo- 
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where 6(. . .) is the Heaviside step function, the zero- 
temperature limit for Fermi-Dirac distribution. The local 
gap equation is 
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The trap introduces the new length scale called 
Thomas- Fermi radius defined as Rtf = \ ■ Note 

further that in the BEC regime, chemical potential fj, is 
already negative at the center of the trap and hence /i(r) 
does not vanish. We also note that in the deep BEC 
regime fj, — —Ef)/2 j23| , the molecular binding en ergy 
with El, — l/2mrO? . Thus we impose the condition |24| . 
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The zero subscript indicates that the quantities are for 
zero polarization. To investigate the system numerically 
we define the dimensionless quantities A(r) — A/Ep, 
fi{r) = fi{r)/Ep, h{r) = h{r)/Ep where we choose Ep = 
(GNy^^hriT with ujt = y^^i + We, also, normalize 

the distance in the trap x = v/Rtfo and define kp hy 
the relation Ep = kp/2rri. Hence 



(11) 
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where we have expressed the binding energy Ei, in Ep 
units as Ei, = Ep / {kpof' However for the system with 
population imbalance we have 



li{x) = p.- Ao + 
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h{x) = h — rj [ fiQ + 



[kpaf 
1 



[kpo] 



(12) 



where dimensionless quantity 77 = [uj^ — uj^^)/{uj^ + w^) 
controls the trap asymmetry of the Fermi gas. 

Now the equation for total number of atoms and pop- 
ulation difference in dimensionless form are given by 



N = j d^xn{ji{x),h{x)), 
AN = / d^xm{fi{x),h{x)). 



(13) 
(14) 



The system for a given coupling strength and polar- 
ization P — AN/N is investigated numerically in the 
following manner: first is calculated by setting P = 0. 
Then using Eq. [12] together with number and popula- 
tion imbalance equation, /x and h are calculated. In the 
BEC limit the order parameter (A) and density for com- 
posite bosons or molecular density are related [2^. By 
calculating the local composite boson density {rim) and 
magnetization (m), the various phases are identified. 
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FIG. 3: (Color online) The three radii Rfi (outer boundary of 
unpolarized superfluid), Rtf (outer boundary of BPl phase) 
and Rf2 (outer boundary of N phase) plotted as a function 
of polarization P for various values of the trap asymmetry 
parameter 77. (a) 77 = —0.9, (b) 77 = —0.5 (c) rj — 0.5, (d) 
rj = 0.9. All the radii are measured in units of Rtfo (outer 
boundary of superfluid unpolarized cloud) . 



RESULTS AND DISCUSSIONS 

We choose experimentally accessible {kpa)~^ = 2.0 
for the interaction strength. The phase at each spatial 
point of the trap is determined by the local chemical po- 
tentials /i(r) and h{r) (see Eqs. [T2|) mapping it to the 
corresponding point in the phase diagram. As radius is 
increased (fJ-ir), h{r)) moves towards left in the phase di- 
agram forming a line segment. 

We find three different phases in the cloud. At the 
center, superfluid core where the population of the two 
components are equal i.e. unpolarized superfluid (BCS 
SF), then an intermediate polarized superfluid (BPl) 
shell where fermion quasiparticles and and composite 
bosons coexist and finally outer rim of majority com- 
ponent. This leads to three radii characterizing the shell 
structure: Rfi where ^ and m{r) — forming 
boundary for BCS SF phase, Rtf above which = 
and Rf2 above which m(r) = 0. 

The three radii for the system without trap asymme- 
try 77 = as a function of polarization P together with 
density profiles showing nm{r) and m(r) are shown in 
Fig. [21 The shell structure consist of BCS SF phase for 
r < Rfi , BPl phase for Rfi < r < Rtf- and finally 
polarized normal (N) state for Rtf < r < Rf2- 

We next consider the system with trap asymmetry 
characterized by the dimensionless parameter ij = (oj^ — 
a;J)/(a;| + wj). The positive (negative) 77 value indicates 
that majority (minority) component is more tightly con- 
fined harmonically than the minority (majority) compo- 
nent. 
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FIG. 4: (Color online) Density profiles at P = 0.65 for differ- 
ent values of the trap asymmetry parameter 77. (a) 77 = —0.9, 
(b) 7? = -0.5 (c) 77 = 0.5, (d) 77 = 0.9. The molecular den- 
sity rim and the magnetization plotted as a function of radius 
measured in units of fc|. and Rtfo respectively. 



The three radii with different trap asymmetry param- 
eter 77 as functions of polarization P are shown in Fig. [S] 
The value 77 ± 0.9 corresponds to the situation when one 
of the component is very strongly confined. We start 
with 77 = —0.9 corresponding to oj^ 4; Wj. The BPl shell 
here is very narrow and overall size of the cloud ( charac- 
terized by Rf2) is much larger than the superfluid cloud 
without population imbalance (the cloud size is measured 
in units of Rtfo)- As we increase 77, the BPl shell grows 
in size, however, size of the cloud decreases. The window 
of polarization for which BPl phase forms the superfluid 
core starting at the center of the trap increases becom- 
ing maximum at 77 = 0.9. The BPl phase forms the core 
beginning at P = 0.76, which should be experimentally 
feasible. We also present the density profiles for same set 
of 77 at P = 0.65 in Fig. HI We note that as 77 increases, 
size of rim representing density of the composite bosons 
shrinks but becomes more dense. It also exhibits the 
large cloud sizes for tightly confined minority as noted 
above. 

All these features can be explained by analyzing 
(/i(r), /i(r)) variations for each value of 77 in the phase 
diagram. To this end, we rcplot the phase diagram en- 
larging the BPl state. The line segments representing 
above mentioned variations are also shown (Fig. [5]). We 
note that 77 with positive (negative) value has a positive 
(negative) slopes with zero value for 77 = 0. 

Note further that BPl to polarized N transition is sec- 
ond and first order for positive and negative set of chosen 
values respectively. These transitions can be detected via 
density profiles in the experiments. 

The 77 = —0.9 line segment traverses small region of 
BPl phase making second order transition to polarized 
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FIG. 5; (Color online) The part of phase diagram with the 
BP-1 region enlarged. The solid (red) line indicates the sec- 
ond order phase transition and dot dashed (black) line shows 
first order transition. Superimposed are the jj. — h variations 
as line segments for various values of the trap asymmetry pa- 
rameter -q at P = 0.65 

N while 77 = 0.9 segment traverses larger region in BP-1 
phase before making first order transition to polarized 
N. The order of transition can be detected in experi- 
ments via spatial discontinuities which vanish for second 
order transition. As 77 is increased from rj — —0.9 line 
segment increasingly have larger portion in BP-1 region. 
This explains why the BPl shell expands in size as 77 is 
increased. Note futher that owing to their negative (pos- 
itive) slopes, rj with negative values have longer (shorter) 
excursion into the polarized N state before encountering 
the vacuum explaining the larger (smaller) size for their 
clouds. Since all the atoms are integrated across the trap 
to conserve the atoms, the corresponding atom density 
distributions are also affected accounting for increased 
(reduced) number of atoms in BPl phase for uj^ > uJi 

CONCLUSION 

In conclusion, we have studied in this paper the 
asymmetrically trapped and population imbalanced two- 
components Fermi gas in the strongly attracting BEC 
limit at T = 0. Using the local density approximation 
(LDA), we calculated the the shell radii of various phases 
in the trap as a function of polarization and trap asym- 
metry. Compared to symmetric trap case (77 = 0), we find 
that when the majority component is tightly confined the 
gapless superfluid shell (BPl) increases in size. The po- 
larization threshold to form the polarized BPl superfluid 
at the core is reduced for a given interaction strength in 
this case. However, when minority are tightly confined 



unpolarized superfluid is favored with BPl phase form- 
ing a narrow shell. We explained these features using the 
phase diagram. 
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